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Abstract. We exploit the 14 gamma-ray bursts (GRBs) with known redshifts z and 

r-t ■ the 7 GRBs for which there are constraints on z to determine the GRB rate -Rgrb(z), 

H~S . using a method based on Bayesian inference. We find that, despite the qualitative 

qq ■ differences between the observed GRB rate and estimates of the SFR in the universe, 

£*\j ' current data are consistent with Rgbs,(z) being proportional to the SFR. 

^H ■ 

> ■ 1 Introduction 

There is increasing evidence that GRBs are due to the collapse of massive stars 

(see, e.g., [1] for a discussion of this evidence). If GRBs arc indeed related to the 

f— «) ■ collapse of massive stars, one expects the GRB rate to be roughly proportional to 

7—1 \ the SFR. However, the observed redshift distribution of GRBs differs noticeably 

^^ ■ from that of the SFR: the observed GRB redshift distribution peaks at z « 1 

(— I \ and few bursts are observed beyond z ~ 1.5, while the SFR peaks at z ss 2 and 

O ,■ 10-40% of stars are thought to form beyond z = 5 (see, e.g., [2,3,]). 

' \ However, observational selection effects play an important role in determining 

5^ ■ the observed redshift distribution of GRBs. The important question is therefore 

"^ \ whether or not the discrepancy between the observed GRB redshift distribution 

& ■ and the redshift dependence of the SFR is entirely due to selection effects; i.e., 

h> \ is the GRB rate roughly proportional to the SFR after taking observational 

•^ ■ selection effects into account? We address this question in this paper. 



2 Method 



We adopt a Bayesian approach. We calculate the likelihood of the data given the 
model, and convert it to a posterior distribution on the model parameters. We 
assume a very general model for the GRB rate, and a power-law model for the 
intrinsic GRB photon luminosity distribution (we assume that the amplitude 
and the power-law index of the photon luminosity distribution does not evolve; 
we relax this assumption in future work). We determine the efficiency with which 
BcppoSAX and the IPN detect GRBs as a function of peak photon flux P by 
comparing the BcppoSAX and IPN peak photon flux distributions to that of 
BATSE. We fit the model jointly to the peak fluxes and redshifts of the 14 
GRBs with known z, and the 7 GRBs for which there are constraints on z. 
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Fig. 1. Comparison of the best-fit models 
to the differential distribution of peak fluxes 
of BATSE, 1PN and BeppoSAX bursts, and 
the cumulative peak flux distributions for 
the three experiments. 
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Fig. 2. Efficiency e(P) with which 
BATSE, 1PN and BeppoSAX detect 
GRBs as given by the best-fit models to 
the differential distribution of peak fluxes 
of BATSE, IPN and BeppoSAX bursts. 



We write the rate of GRBs that occur per unit redshift and luminosity as 



dN/dzdL N = p{z)f(L N ), 



where 



p(z) = Rgkb(z;P,Q) x (1 + z)- 1 x ^r(zf(dr/dz) 
is the rate of GRBs that occur at redshift z, 

p 



Rgkb(z;P,Q) 



t(z) 



t(0) 



exp 



1 



*(£) 
t(0) 



(1) 
(2) 

(3) 



is the rate of GRBs that occur at redshift z per unit comoving volume (see Q ) , 
t(z) is the elapsed time since the Big Bang, and 



■-P 



f(L N ) = L~ p x 0{L N - L min )0(L„ 



■>N) 



(4) 



is the intrinsic photon luminosity distribution of GRBs. Thus the model has five 
parameters: P, Q, /3, L min , and L max . 

We write the efficiency with which GRBs with known redshifts are found as 
e(z, P) — € z (z) x esTtp(P), where e z (z) is the efficiency with which the redshifts 
of GRBs are determined from optical observations once they are detected by 
a 7-ray burst instrument. We take e z (z) = 1 for GRBs whose redshifts were 
determined by detection of an absorption-line system in the optical afterglow 
of the burst and e z (z) = 0(1 — z) + Q(z — 2.5) for GRBs whose redshifts were 
determined by measuring emission lines in the spectra of the host galaxy. The 
latter expression accounts qualitatively for the difficulty in measuring redshifts 
when the H a and 0[ff] emission lines from host galaxies do not lie in the visible 
spectrum. The quantity cst is the "stereo-temporal" efficiency that accounts for 
limitations of exposure in time and solid angle and ep(P) is the efficiency with 
which BATSE, the IPN and BeppoSAX detect GRBs as a function of peak flux 
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Fig. 3. Probability contours in the 
(P,Q)-plane for the GRB rate param- 
eters (P, Q) found from fitting to the 14 
GRBs with known z and the 7 GRBs 
with constraints on z. The solid curves 
correspond to the 68% and 95% prob- 
ability contours. Also shown are the 
(P,Q)-values corresponding to no space 
density evolution (+), the Madau et al. 
SFR [4], and the Rowan- Robinson phe- 
2(i - nomenological model fit to IR, optical 

5 10 and UV data [5]. 

P 

P. Figure 1 compare the best-fit models of ep and the cumulative peak flux 
distributions of BATSE, the IPN, and BeppoSAX, respectively. Figure 2 shows 
the best-fit ep for each of the three experiments. 
The likelihood function is then given by 

N 



where 



C = cxpl - / dzdP /j,(z,P)e(z,P) > JJ/i(^»,-P») , 

^ J i—l 



(5) 



L 



,( Z ,P) = J o dL NP ( Z )f(L N )S^- 47rr{z ^ i + z)a 

= p(z) x / (47rr(z) 2 (l + z) a P) x 4nr(z) 2 (l + z) a (6) 

is the expected number of events observed within dz dP of (z, P). The quantity 
a is the burst spectral index, which we set equal to one in this work. By an ap- 
plication of Bayes' Theorem, we now regard C as an (unnormalizcd) probability 
distribution on the model parameters. 

3 Results 

Figure 3 shows 68% and 95% probability contours for the GRB rate parameters 
P and Q. Also shown on the plots are the best-fit SFR models of Madau et 
al. [4] and Rowan- Robinson [5]. The SFR models lie at about a 68% excursion 
from the best-fit GRB rate model. Thus we find that, despite the qualitative 
differences that exist between the observed GRB rate and estimates of the SFR 
in the universe, current data are consistent with the actual GRB rate being 
approximately proportional to the SFR. 
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